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Abstract 

In this paper we present a simple analytic expression for the (spin- 
averaged) neutrino-nucleon cross section for ultra-high energies at 
twist-2, obtained as the asymptotic limit of our previous findings. This 
expression gives values for the cross section in remarkable numerical 
agreement with the previous numerical evaluation in the energy region 
relevant for forthcoming neutrino experiments. Moreover, we discuss 
the role and the relevance of saturation and recombination effects in 
our approach, in comparison with other recent suggestions. 
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1 Introduction 


Neutrino astronomy holds enormous scientific potential and prospects for 
its development are much better at high energies since the neutrino-nucleon 
cross section and the angular resolution increase with energy [l|. These fea¬ 
tures give the opportunity to use large natural target media like ice and 
atmosphere as detectors. The ultra-high energy neutrino-nucleon cross sec¬ 
tion will become soon an important ingredient in the interpretation of the 
results in these experiments. 

Above the energy at which the neutrino interaction length is approx¬ 
imately equal to the diameter of the Earth, ~ 40 TeV, no experimental 
constraints exist on the neutrino cross section. However, parton distribution 
functions fitted to the HERA data and Standard Model constrain theoretical 
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obtained with different models and all cross sections 


predictions 

are remarkably consistent at the highest energies [^. In some models 
the introduction of non-linear screening effects produces only mild changes 
and the aforesaid remarkable agreement between cross sections survives; ac- 
cording to other models [3, all-twist formulation of QCD evolution equa¬ 
tions [ill, [l2, ll3| entails a drastic change on the cross sections in the region 
where unitarity effects become important. Geometric scaling that is a 
consequence of the all-twist formulation of QCD, and a precise form of the 
dipole cross section 1^ in the geometric scaling region are assumed in the 
approach of Ref. |]i|. We will try to explain roughly these concepts since 
they will be important in the following. 


In the color dipole picture m the virtual photon, or the gauge bosons 
W and Z for neutrino scattering, creates a qq pair, or a “color dipole”. At 
high energies, or at small x, the exchange of gluons between the nucleon and 
the color dipole becomes more and more important and the gluon density 
in the nucleon increases with energy. The quark-antiquark dipole has zero 
color charge and the interaction with the gluons in the nucleon will depend 
on the dipole size. If the size is very small the dipole will not interact with 
gluons, but there will be a length Ro{x) = 1/Qs{x) such that, for dipole sizes 
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greater than -Ro(3;), the scattering cross section will be perceptible. i?o(a^) 
is called the saturation radius and Qs{x) the saturation scale. Geometrical 
scaling refers to the dependence of the dipole cross section a{x,r), where 
r is the transverse separation of the quarks in the qq pair, from only one 
dimensionless variable t/RqIx). As a consequence the p cross section, for 
example, becomes a function of one dimensionless variable r = Q‘^R^{x) (l4| . 
At small X, and hence at high gluon density in the nucleon, Ro{x) is small 
and the saturation scale is large. 

We return now to the neutrino-nucleon cross section. According to Ref. 
ful . at neutrino energies > 0(10^^) GeV geometric scaling holds and 
neutrino-nucleon cross sections are enhanced by a large factor. A higher cross 
section in this energy region could have important consequences for neutrino 
astronomy Q- Hence, a comparison with approaches, where geometric scal¬ 
ing and saturation scale have a different validity range and interpretation, 
becomes interesting. 

In a series of papers jl8|, ll9| modihed evolution equations were suggested 
that include twist-4 gluon recombination corrections. These equations are 
quoted as “modihed DGLAP equations” because of their similarity with the 


original evolution equations j2l 


Refs. 


00 


10 : 


gluon recombination is evaluated in 


with the same technique adopted in Ref. 


211 at twist-2. The 


physical picture of the deep inelastic scattering (DIS) process is different from 
the color dipole one, since now the virtual gauge boson explores the parton 
distribution in the nucleon. In the leading logarithmic approximation, and at 
twist-2 level, the t wo p ictures give the same results but, at higher twist, this 


equivalence is lost since approximations are different in different frames. 

22| . the most important reason for this difference is that 


According to Ref. 

the color dipole approach extracts the splitting probabilities incoherently and 
neglects the coherence amon g di fferent subpartonic amplitudes. 


It has been shown in Ref. 


23| that the Balitsky and Kovchegov non-linear 


evolution equation [11|, |2^ leads to saturation of the scattering amplitude, 
but does not necessarily unitarize the total cross section. The violation of 
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unitarity depends on the nature of the target and the reasons for its ap¬ 
pearance can be seen both in the target rest frame or from the point of 
view of the evolution of target helds. Interactions between the dipoles in 
the projectile wave-function are neglected in the color dipole approach and, 
on the other side, the evolution in the target is driven from incoherent color 
sources. Unitarity is violated since the dipole interacts with the long range 
Coulomb held created by a large number of incoherent color sources in the 
target. The proof in Ref. Qi holds for a strong interacting particle colliding 
on a hadronic target, but for DIS the situation changes for the worse since 
the DIS cross section acquires an extra power of the rapidity. As noticed in 


Ref. j^, the condition that the total color charge must be zero, in a region 
of hnite size (e.g. the proton), introduces correlations, and coherence, among 
the sources of the color charges and can lead to a unitary evolution. These 
considerations justify the interest for a comparison between the two different 
approaches to saturation and its consequences for neutrino astronomy. 


In our previous paper [9| we have obtained a simplihed solution of the 
non-linear evolution equations of Refs. EQ at small x. In this paper we 
will hrst show that it is possible to simplify further the integrals leading to 
the charged current neutrino-nucleon cross section (Section 2). The answer 
at twist-2 will be analytical and take a very simple form. In this simpli¬ 
hed approach it becomes easy to prove the well known statement that, at 
asymptotic energies, only the values of 


tion P,l2^ and an estimate of the neglected terms will be given (Section 3). 

With the parameters, in the input parton distributions, hxed at the values 
obtained in Ref. [£ L we will then discuss the saturation phenomenon in the 
approach of Refs. 18|,|l9|. Similarities and differences with the approach in 
Ref. El will hnally emphasized (Section 4). 


contribute to the cross sec- 
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2 Asymptotic form of the neutrino-nucleon 
cross section 


The starting point is the inclusive, spin-averaged cross section for the neu¬ 
trino interaction with an isoscalar nucleon target, N = (neutron-|-proton)/2, 
in the process 


iy^{k) + N{p) ^ H{k') + X{p') , (1) 

where parentheses enclose the four-momenta of the particles participating to 
the scattering. The transformation of the neutrino to a charged muon labels 
the event as a “charged current” event and the charged current cross 
section can be expressed in terms of the nucleon structure function as 

( da y ^ GlME, / y 

\dxdy) TT + 


X 






.( 2 ) 


For anti-neutrino charged-current processes one must change the sign in 
front of while the changes necessary in order to describe neutral-current 


0. In Eq. Q = 


neutrino interactions can be found in the literature 
Fa x,Q‘^) for i = 1, 2, 3, Gp is the Fermi constant, M is the nucleon mass 
and Mw is the hF-boson mass. The scaling variables x and y are defined as 

,2 


X = 


T 


2p ■ q 


y = 


p-q 

p ■ k 


( 3 ) 


and = —q^. 

At high energies the relation between the variable y and the Bjorken 
variable x can be approximated as 


^ a:(s — M^) 




( 4 ) 


where s = {k + p)'^ is the square of the c.m. energy for the neutrino-nucleon 
scattering. The laboratory neutrino energy Ey = [s — M^)/{2M) is approx¬ 
imately equal to s/{2M) in the region we consider. 
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The main purpose of this paper is to evaluate the asymptotic behavior of 
the total cross section for the charged-current neutrino-nucleon process 


a 


= dx dy 


r da y 


( 5 ) 


\dxdy J 

In the following, we will limit ourselves to the leading order corrections to the 
simple parton model and hence all parton model formulas remain unchanged 
except that the parton distributions depend now on x and and not only 
on X. In particular, the Callan-Gross relation, F 2 = 2xFi or Fl = 0, holds in 
leading order. By imposing a lower cut in = Qq, we rewrite Eq. 

as 


a 


vN {CC) 


1 


2ME, 


dQ^ 


dx f da \ 


JQl 




X 


\dxdy) 


where y = Q^/(xs) in the differential cross section, according to Eq. (jH). 
Then, the F 2 contribution to the total cross section can be written in the 
form 




a^N + 


/*5 


2n 








dx 1 -I- (1 — (xs))^ 


'Q^/s 


X 


■F,yx,Q^), ( 6 ) 


where we have neglected the term —Mxy/{2E) in Eq. Q. Since we are 
mainly interested in the asymptotic s behavior, some simplifications are pos¬ 
sible and, when s is much larger than all the scales appearing in Eq. ©, in 
particular s S> 

1. the contribution of xF^ can be neglected and ~ a'^^ ~ a^^] 

2. the inequality s ^ holds because the factor -|- M^)Y 

limits the integration region: as will be shown later, the upper limit 
of the integral becomes proportional to M^. 


As in Ref 


S. 


we write the isoscalar structure function in terms of the 


parton distribution functions and, at leading order, we get 

F 2 {x) ~ X u{x) + X u{x) + X d{x) + X d{x) -|- 2x s{x) + 2x c{x) + 
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where the dots stand for the b- and t-quark PDFs and we have assumed 
s{x) = s{x) and c{x) = c{x). We denote, in the following, by fq{x,Q^) the 
sea quark distribution xS{x, at twist-2 and by fg(x, Q^) the gluon distri¬ 
bution xG{x, Q^) in the same approximation and use the notation Q^) 

for the sea quark distribution modihed by the introduction of gluon recom¬ 
bination at twist-4 [l^ . 

Setting u = Q'^js, we can write the last integral in Eq. in the form 


dx 

, X 


^ u u 

X 2x‘^ 




(7) 


which is a Mellin convolution that, in moment space, becomes 

dx 

X 

1 




M 


n — 1 


j + 0(n-l)^ //"“(n.Q"), 


( 8 ) 


since, as noticed before, large contributions are strongly suppressed by 
the factor (1 -|- and u can be considered small. The result in 

Eq. (jHI) can be obtained by considering the relation 

dx 


u^-^du / —Ml f-) M^ix) = Mi(n)M 2 (n) , 
.L, X \xJ 


where 


Mi(n) = 


'0 


dte-^ 
1 1 


1 — t -f- 


n 


n 


2{n -|- 1) n — 1 4 


^ -^ + 0(n-l) 


and 


M 2 (n)= / 


With the dehnition 


n 




9{u,Q^) , 


(9) 
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Eq. (jni) becomes 

3 Twist-2 contributions to the cross section 


To begin with we consider our approximate twist-2 solution for the structure 
function F 2 {x, Q^), that is the hrst term fg in 

//“" = f, + T,, 


where Tq represents the twist-4 gluon recombination corrections. It is impor¬ 
tant to ensure the accuracy of the approximation oa in the simplest case 
and to e^lore the possibility of further simplihcations. From our previous 



work 
also in Ref. 


, where the method introduced in Refs, 
was adopted, we have 
1 


27, 


281 129 I. Iso I and used 


fq{n,Q ) = 


n 






-dggt 


k+l 




fc =0 


(fc -|- 1)! \n — 1 


( 11 ) 


where 


f = In 


O^siQl) 


dgg = -121 po, d_{l) = 16//(27/3o) and d+(l) = 1-0 20//(27/3o), with po = 
11 — 2//3 and / the number of flavors. Introducing the new variables 

CTqi 


(Jn, - 2 


-dggt\n.{l/u), Pu = 


21n(l/M)’ 

we hnd the following expression for the twist-2 contribution to g{u, Q^): 


( 12 ) 


g^‘^\u,Q^) ~ A- 


u 4 




+ A+ 


do(^u) ^Pudl(o'u) 




(13) 
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The simplest contribution to Eq. m coming from Q'^) is 


where we can put exp[—= [ln(QQ/A^)/ln(Q^/A^)]'^-(^). A basic inte¬ 
gral appearing in Eq. m is then the following: 


hid) 




1 

[ln(QVA2)]<^ • 


The proof that, when s —> cx), 


(14) 


hid) ^ M, 


w 




( 15 ) 


neglecting terms proportional to [ln(M^/A^)]“^“'^, is rather long, but im¬ 
portant and will be presented in Appendix A. 

It is not difficult to generalize the proof of Eq. m to an expression of 
the form 

“ IqS ^ • 

where j[z] is a function that can be expanded in powers of 1/z, and prove 
that, at twist-2, 


a 


vN ( 2 ) 


27r 


A- In 




- 4 '" 


—d— (l)i 




( 16 ) 


where 
£ = In 


O^siQo) ^ 

Ws{M^w) 


a = 2 


-dggi\n{s/M‘^), p = 


a 


2\n{s/M^) 


We notice that the expression for the twist-2 cross section is explicit. In 
Table 1 and Fig. 1 this approximation is compared with the numerical de¬ 
termination obtained in Ref. P| in the case of absence of recombination (see 
















Fig. 4 of that paper). The values used for the parameters Ag and Ag are 
1.040(36) and 0.548(28), respectively The approximate cross section given 
in Eq. (fT^ nicely matches the numerical determinations of our previous 
work and those of Refs. Bill- 

When the gluon recombination term is present, any attempt to simplify 
the problem becomes much more intricate but, as we will see in the next 
Section, the discussion of the saturation limit can be done on the basis of a 
simpler approach. 


s [GeV2] 

a^N [cm2], ^ 

auN [cm2], 

10^ 

9.75(39) X 10-35 

1.027(38) X 10-34 

10^ 

4.13(17) X 10-34 

3.65(15) X 10-34 

10^ 

1.336(57) X 10-33 

1.112(46) X 10-33 

10^ 

3.75(16) X 10-33 

3.20(14) X 10-33 

10^ 

9.67(42) X 10-33 

8.61(37) X 10-33 

IQio 

2.33(10) X 10-32 

2.164(94) X 10-32 

IQii 

5.34(23) X 10-32 

5.11(22) X 10-32 

1012 

1.167(51) X 10-31 

1.143(50) X 10-31 

IQi^ 

2.45(11) X 10-31 

2.44(11) X 10-31 

IQi^ 

4.97(22) X 10-31 

5.02(22) X 10-31 

1015 

9.80(43) X 10-31 

9.96(43) X 10-31 


Table 1: Comparison between the twist-2 contribution to the neutrino- 
nucleon cross section according to the numerical results obtained in Ref. B 
(2nd column) and the asymptotic approximation given in Eq. ()16|1 (3rd col¬ 
umn). The errors come, in both cases, from the uncertainties of the A^ 
parameters. 


^We remind that 


f{Ag + AAg/9)/9 and A = Ag, for flat initial conditions 
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4 Saturation and recombination scales 


The saturation scale Q|(x) indicates the saturation limit and is usually de- 


hned as 


32 


d xG{x, 


d lng2 

, equivalently [l^ . 


= 0 and 


d xS{x, 


Qi 


d lng2 


= 0 




1^5 = 


non-linear terms 


linear terms 


= 1 , 


which means that the non-linear recombination effect in the MD-DGLAP 
equation fully balances the linear splitting effect. The recombination scale 


g|., introduced in Ref 


0 , 


is defined as 


non-linear terms 


linear terms 


= as[QUx)] , 


Q? 


which means that, near this scale, the higher-order recombination contri¬ 
butions cannot be neglected and should be included in the evolution, thus 
making the evolution of the parton distributions from to g| much more 
complicated. 

Since, according to Ref [l^, saturation and recombination appear at very 
low X we are justihed to use approximate relations like 


dF^^{x,Q^) _ df^^'‘\x,Q'^) dfg{x,Q^) al 


17 


+ (17) 


dlng2 dlng2 (ilng2 ' g 2 " 32''^ 

(see, for example, Eq. (64) in Ref and Appendix B). The saturation scale 
will be consequently dehned from the equation 


dfq{x,Q'^) 
d In 


IIJL f2(j. /n2'i 
32 Jg I"'') V ) 

while the recombination scale satishes 

dfq(x,Q^) 

dlnQ^ 




(18) 


17 K 




— Qio 


(19) 
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We can use in these equations the values of the parameters obtained in the 
£t of ZEUS PDF in order to have an idea of the behavior of the scale 
and more importantly of the scale with x. This is an interesting point 
since twist-4 recombination formulas hold near the recombination scale but, if 
we approach the saturation scale, higher order recombination contributions 
become signihcant In other words, we identify the region where our 

formulas can be trusted. From Eqs. (HHD and (P|) . we can build a numerical 
table (see Table 2 and Fig. 2), using the results of our previous work [^, 
where, in particular, the K parameter was set to 0.013. 


Q2 [GeV^] 

Xf'ec 

^sat 

5 

1.35 X 10-® 

3.27 X 10-1^ 

10 

8.43 X 10-9 

1.79 X 10-1^ 

50 

7.49 X 10-1° 

9.14 X 10-13 

100 

1.96 X 10-^0 

1.82 X 10-13 

200 

4.59 X 10-^1 

3.20 X 10-1^ 

500 

5.87 X 10-^2 

2.78 X 10-13 

1000 

1.14 X 10-^2 

4.00 X 10-13 

2000 

2.06 X 10-^^ 

5.39 X 10-1^ 

3000 

7.37 X 10-^^ 

1.62 X 10-1^ 

4000 

3.52 X 10-^^ 

6.82 X 10-13 

5000 

1.97 X 10-^^ 

3.46 X 10-13 


9.98 X 10-^® 

1.57 X 10-13 


Table 2: Recombination and saturation scales, dehned according to Eqs. (HHD 
and (P|) . determined using the PDFs of our previous work 

From this table we realize that 

1. the function is quite similar to the one obtained in Ref. with 

different values of the parameters (in particular the value of K is quite 
different in our approach); 
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2 . the recombination scale agrees with our hndings for the slope of 

3. the evaluation we did of the neutrino-nucleon cross section is safe since 

the small-x limit considered is well above the a:-value associated with 
the recombination scale at ~ in other words higher order 

twists, besides twist-4, are not important. 


This last point gives sense to a comparison of our results with those 


io|, where all twists were resummed, in the region of energy we 


of Ref. 
considered. 

As far as anti-shadowing effects are concerned, they are important at val¬ 
ues of much smaller that , therefore they are completely negligible 

in our analysis. 

In conclusion we hnd that our approach is sound and reliable. It has 
many points in common with the analysis of Ref. saturation presents 
itself at very small x and is not detectable, antishadowing is present in both 
approaches, but at different values of Q^. Enhancement of ultra-high energy 
neutrino-nucleon cross section is not required in our calculation and this is 
due to the correlation among color sources in the target. Such correlation, 
and coherence, is absent in the color dipole picture. 
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s [GeV"] 


Figure 1: Comparison between the twist-2 contribution to the neutrino-nucleon 
cross section according to numerical result shown in Fig. 4 of Ref. (squares) and 
the asymptotic approximation given in Eq. (Illil) (continuous lines, representing the 
upper and the lower bounds at la level). The uncertainty comes, in both cases, 
from that of the parameters. 
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[GeV ] 



Figure 2: Recombination (left) and saturation (right) scales, defined according to 
Eqs. m and (HU . determined using the PDFs of our previous work |^. 
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A Appendix: Proof of Eq. (1 151} 


With the change of variable t = /h?) we can rewrite the integral (d 

in the following forms 

r-ln(s/A2) 


Ti{d) = 




t dt 


iin(Qg/A^) (e* + M^/A^f 

f-d 


A2 

A2 


T\yf2 f\n{s/AA 


t dt 


/ln(Q2/A2) e* + Af^/A2 A2 iin(Q 2 /^ 2 ) (e* + M^/A2) 
Hn(./A2) dt _ ^ dt 

/ln{Q2/A2) 1 + M^e“VA^ A2 Jin(Q2/A2) (1 + 


Setting z = M^/A^, we have 


Ji(h) = M^z{l + z—j 


d \ t-^e-^ dt 


'ln(Q2/A2) 


1 + ^e“* 


d 


= M^z [1 + z—] C{d,z) . 


C{d, z) differs from the integral 


hn(Qg/A2) 


dz 


t ‘^e ^ dt 

~ 


(A.l) 


1 + ze 


by terms vanishing faster than 1/s in the asymptotic region for the variable 
s. This resnlt can be easily obtained by expanding in series the integral in 
Eq. (EH with respect to its npper limit. Moreover, since z is a large nnmber 
(In 2 ; ~ 12 if A = 0.19 GeV) another approximation becomes possible and, 
neglecting terms proportional to [ln(^)]“^ with respect to a constant term, 
we have 

The integral (Ha can be expressed as an inhnite snm 

00 

c{d,z) ~ 

n=o do 
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= 4^(1(A-3) 

n=l 

and finally 

Cc^--r{l-d)F{-z,l-d) , (A.4) 

z 

where F{z,2) is the Euler’s dilogarithm. The analytical continuation of the 
series 

F(-z,l-d) = J2 

n=l 

is given by the Joncquiere’s relation S that, in our case, becomes 

F{-z, l-d) = 1 - d) + (^d, j , 

The variable z = M^/A^ is very large and an asymptotic expansion for C 
follows from the asymptotic expansion of the generalized Zeta function for 
z ^ oo, 


Since 


and 


CM, 


ln(—z) 
27ri 


T{d) 


r(d- 1 ) 


ln(—z) 
271 i 


l-d 


-d 


+o(i>"Wr‘T 


d 


li{d) = d)z\^ + z—j^{-z,l- d,l) 

= -Ml,T{l-d)F{-z,-d) 


-F{-z,-d) = -e-^^‘^/^r{d){-i\nz)-‘^sm{77d) , 

TT 

Eq. (fTH|l follows at once. 


(A.5) 


(A.6) 

(A.7) 
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B Appendix: Derivation of Eq. (1 171} 

According to the modified DGLAP equations j^, we have 


Q^) 

dlng2 


dfgjx, Q^) ^ 

dlng2 +Q2 


.1/2 


y 


F, 


99 


X 


y 


y 

(//""(!/.(?/) 




Q^) , 

(iing2 g2 


r'tc 


dy 


X 


99, (B.l) 

L^x/2 y \y. 


-1/2 


y \yj 


“ dlng2 

Then, using Eqs. (28), (29) and (32) of Ref. and recalling that we put 
Ki = K 2 , we easily get 


Rq{x,Q‘^) 


M 


n^l 


M- 


d/(") - q/(")l h(n) - FS{n) h(n) 


<19 

131 


99 


1813 \ , , , 17 , , , 


(B,2) 

(B.3) 
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